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Abstract 

Given a planar graph G on n vertices and an integer parameter r < n. an r division of G with few 
holes is a decomposition of G into 0{n/r) regions of size at most r such that eacli region contains at most 
a constant number of faces that are not faces of G (also called holes), and such that, for each region, the 
total number of vertices on these faces is 0{^/r). 

We provide a linear-time algorithm for computing r-divisions with few holes. In fact, our algorithm 
computes a structure, called decomposition tree, which represents a recursive decomposition of G that 
includes r-divisions for essentially all values of r. In particular, given an increasing sequence r = 
{ri,r2, ■■■), our algorithm can produce a recurstve r-division with few holes in linear time. 

r-divisions with few holes have been used in efficient algorithms to compute shortest paths, minimum 
cuts, and maximum flows. Our linear-time algorithm improves upon the decomposition algorithm used 
in the state-of-the-art algorithm for minimum st cut (Italiano, Nussbaum, Sankowski, and Wulff-Nilsen, 
STOC 2011), removing one of the bottlenecks in the overall running time of their algorithm (analogously 
for minimum cut in planar and bounded-genus graphs). 
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1 Introduction 



Separators decompose a graph in a balanced way into two subgraphs with a hmited number of vertices in 
common. Separators are often used in efficient algorithms using a divide-and-conquer strategy LT80, Chu9Ql 
IShm96|. IRHOOj . Graphs with smaU recursive separators include planar graphs [UngSl, iLT79.. .Dji82, ,Mil86| 
IGM901 [ST961 IDV97j, bounded-genus graphs [Dji85| IGHT841 IKel06j . mino r-free graphs |And86l IAST90I 
IPRS941 IRW091 IBLRlOl iKRTOl IWNllj . and graphs with bounded tree- width [HilTel IREMI RP89] . Further- 
more, for graphs of all these classes, separators can be found efficiently, often in linear time. For planar 
graphs, experimental results demonstrate that separator algorithms are practical |ADGMd6l lHSW"'"09 . 

Perhaps the most influential result of this kind is the linear-time algorithm of Lipton and Tarjan jLT79] 
for finding a separator of size 0{^/n) in a planar graph with n vertices. Consider the result of using 
this algorithm recursively until each resulting graph has size at most some specified limit r. It is easy 
to show that 0{n/r) subgraphs result, and that the average number of boundary vertices per subgraph 
is 0{y/r). Frederickson jFre87| showed that, with additional care, one can ensure that each subgraph has 
0(-y/r) boundary vertices; he named such a decomposition an r-division, and he referred to the subgraphs 
as regions. The running time of Frederickson's algorithm is 0(n log r -I- (nj yjr) logn). 

Decompositions of such kind have been used in many efficient algorithms for planar graphs, e.g. for 
computing shortest paths |Fre87[ |HKRS97| . maximum flow |JV82j . polygon triangulation |Goo95| . and I/O- 
efficient algorithms |MZ08| . 

For some of these applications, Frederickson's algorithm for r-divisions is too slow. Henzinger, Klein, 
Rao, and Subramanian |HKRS97] . whose linear-time shortest-path algorithm required recursive application 
of an algorithm for r-divisions (with roughly log* n levels of recursion) addressed this by showing that a 
linear-time 0{\pn) separator algorithm could be used to obtain a sublinear-time separator algorithm with 
a worse (but still sublinear) boundary-size guarantee. In a major step forward, Goodrich |Goo95| gave a 
linear-time algorithm that achieves 0{^\fr) boundary size for planar graphs (see |ABT04) for the I/O model). 

However, for some of the more involved algorithms working with planar embedded graphs, it is essential 
that the regions of the division are topologically "nice" in that the boundary of each region consists of a 
constant number of faces, also called holes. Such a division can be found by using small cycle separators 
(Miller ^Mil86j ) instead of just small separators, and incorporating iterations in which the graph is separated 
according to the number of holes. 

Such r-divisions with a constant number of holes were first used in algorithms of Klein and Subrama- 
nian [KS981 [SiiMSj , and subs equently in many other algorithms [FRMI ICRToI ICabl2l IMWNIOI IINSWNIII 
iLSm IKKSTTI IEFN121 |MST2] Q 

Up to now, the fastest known algorithm computing an r-division with a constant number of holes per 
region runs in time 0(nlogr -|- {n/^/r) logn) |INSWNTT| . This makes it one of the time bottlenecks in the 
state-of-the-art algorithms for minimum si-cut and maximum st-flow [INSWNllj and minimum cut |LSllj in 
undirected planar graphs and bounded-genus graphs IEFN12j . Whether such an r-division can be computed 
in linear time was an open problem until the current work. For example, Cabello |Cabl2j remarks that "it is 
unclear if the algorithm of Goodrich \Goo95f can also be modified to use the cycle-separator, and thus obtain 
a linear-time construction of r -decompositions with a few holes." 

Contributions. 

We provide a linear-time algorithm for computing an r-divisions with few holes for any triangulated bicon- 
nected plane graph G and any r (Theorem [ij . 

In fact, the algorithm produces a decomposition tree of G (Theorem [2]), which is a tree that naturally 
represents a recursive decomposition of G by cycle separators, and from which one can read off, in linear 
time, a recursive r-division with few holes, for any increasing sequence r = ri,r2, . . . (Theorem |3]). 

Our linear-time algorithm improves upon the 0(rilogr+ (n/y^) logn)-time algorithm of Italiano, Nuss- 
baum, Sankowski, and Wulff-Nilsen [INSWNll , removing one of the time bottlenecks of the state-of-the-art 
algorithms for minimum s<-cut and maximum st-flow jlNSWNllj . as wefl as minimum cut (LSTT] IEFNT2] . 

^Cabello |Cabl2| requires only that the average number of holes per region be small. 
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Techniques. 



The overall approach of our algorithm builds on that of Goodrich |Goo95j . which is based on Lipton-Tarjan 
vertex separators |LT79I . However, our approach must handle the additional complexity of finding cycle 
separators (see Miller Mil86,), which involves maintaining spanning trees in both the input graph and its 
planar dual, and of bounding the number of holes. We manage both primal and dual trees simultaneously 
by using a primal spanning tree that follows a dual breadth-first-search tree (Section 4.2). The levels of the 
dual BFS tree define connected components, which we maintain in a Component Tree (Section 4.1 1. This 
component tree internally captures the structural connectivity of dual BFS components. 



Outline. 

Precise definitions of r-divisions, holes, recursive divisions, and decomposition trees, as well as formal state- 
ments of our results are in Section |2.5| and Section |2.6| We give a high-level description of the r-division 
algorithm and prove its correctness in Section [Sj 

A key ingredient in our algorithm is a cycle separator algorithm (Section |4]) that can be implemented 



efficiently by using auxiliary data structures (based on dynamic trees and Euler-tour trees, see Section 2.8 1 
In a straightforward computation of an r-division, the cycle separator at each recursive level can be found 
in linear time and the graph can be cut in linear time into two regions along the cycle. In order to achieve 
overall linear time, we aim at suhlinear time bounds for computing the cycle separator and cutting the 
graph along it. In our algorithm, we initialize the data structures once, and then, in subsequent steps of 
the recursion, we reuse what has been computed already (Section [5|, and we call an efficient dynamic-tree 
implementation of our cycle separator algorithm (Section |4.5[) . 



2 Preliminaries 
2.1 Graph Notation 

Let G — (V, E) be a simple graph. For a subset V' of V , we denote by G\V'] the subgraph of G induced by 
V\ i.e., the graph obtained from G by deleting all vertices not in V' and all their adjacent edges. 

For a subset V' of V , we denote by 5g{V') the set of edges uv of G such that v ^V' and u ^ V' . We 
refer to 5g{V') as a cut in G. If V' and V — V both induce connected subgraphs of G, we say it is a simple 
cut (also known as a bond). In this case, V and V are the vertex-sets of the two connected components 
of Sg{V'), which shows that the edges of Sg{V') uniquely determine the bipartition {V , V — V'}. 

A graph is called biconnected iff any pair of vertices is connected by at least two vertex-disoint paths. 

For a spanning tree T of G and an edge e of G not in T, the fundamental cycle of e with respect to T in 
G is the simple cycle consisting of e and the unique simple path in T between the cndpoints of e. 



2.2 Embeddings and Planar Graphs 



We review the definitions of combinatorial embeddings. For elaboration, see also http://plaiiarity.org 
Let £' be a finite set, the edge-set. We define the corresponding set of darts tohe E x {±1}. For e G E, 

the darts of e are (e, -1-1) and (e, — 1). We think of the darts of e as oriented versions of e (one for each 

orientation). We define the involution rev(-) by rev((e,i)) — (e, — i). That is, rev(d) is the dart with the 

same edge but the opposite orientation. 

A graph on E is defined to be a pair {V, E) where ^ is a partition of the darts. Thus each element of 

y is a nonempty subset of darts. We refer to the elements of V as vertices. The endpoints of an edge e are 

the subsets v,v' containing the darts of e. The head of a dart d is the subset v gV containing rf, and 

its tail is the head of rev(c?). 

An embedding of {V^E) is a permutation tt of the darts such that V is the set of orbits of tt. For each 

orbit V, the restriction of tt to that orbit is a permutation cycle. The permutation cycle for v specifies how 

the darts with head v are arranged around v in the embedding (in, say, counterclockwise order). We refer 

to the pair (tt, E) as an embedded graph. 
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Let TT* denote the permutation revo n, where o denotes functional composition. Then (vr*, E) is another 
embedded graph, the dual of (ir^E). (In this context, we refer to {j^,E) as the primal.) 

The faces of (tt, E) are defined to be the vertices of {tt*,E). Since rev o (rev on) = tt, the dual of the 
dual of (tt, E) is (tt, E). Therefore the faces of (tt*, E) are the vertices of (tt, E). 

We define an embedded graph (tt, E) to be planar if n — m + = 2k, where n is the number of vertices, 
m is the number of edges, (j) is the number of faces, and k is the number of connected components. Since 
taking the dual swaps vertices and faces and preserves the number of connected components, the dual of a 
planar embedded graph is also planar. 

Note that, according to our notation, we can use e to refer to an edge in the primal or the dual. 

2.3 Some Properties of Planar Graphs 

Fact 1 (Simple-cycle/simple-cut duality jWhi32] ). A set of edges forms a simple cycle in a planar embedded 
graph G iff it forms a simple cut in the dual G* . 

Since a simple cut in a graph uniquely determines a bipartition of the vertices of the graph, a simple 
cycle in a planar embedded graph G uniquely determines a bipartition of the faces. 

Definition 1 (Encloses). Let C be a simple cycle in a connected planar embedded graph G. Then the edges 
of C form a simple cut 5g* (S) for some set S of vertices of G* , i.e. faces of G. Thus G uniquely determines 
a bipartition {Fq, Fi\ of the faces of G. Let foo, f be faces of G. We say C encloses / with respect to f^o 
if exactly one of /, foo is in S. For a vertex/edge x, we say G encloses x (with respect to foo) if it encloses 
some face incident to x (encloses strictly if in addition x is not part of G ). 

Fact 2 ( |vS47) ) . For any spanning tree T of G, the set of edges of G not in T form a spanning tree of G* . 

For a spanning tree T of G, we typically use T* to denote the spanning tree of G* consisting of the edges 
not in T. We often refer to T* as the cotree of T |Epp03| . 

2.4 Separators 

Definition 2. For an assignment W{-) of nonnegative weight to faces, edges, and vertices of G, we say a 
simple cycle G is a balanced separator if the total weight of faces, edges, and vertices strictly enclosed by C 
and the total weight not enclosed are each at most 3/ A of the total weight. 

(Traditionally balance involves a bound of 2/3. We use 3/4 because it simplifies the presentation.) 

One can reduce the case of face/edge/ vertex weight to the case of face weight. For each vertex or 
edge, remove its weight and add it to an incident face. A cycle separator that is balanced with respect to 
the resulting face- weight assignment is balanced with respect to the original weight assignment. We may 
therefore assume in cycle-separator algorithms that only the faces have weight. 

Let G be a planar embedded graph with face weights. Suppose that no face has more than 1/4 of the 
total weight. Lipton and Tarjan |LT79j show that, if G is triangulated (every face has size at most 3) then for 
any spanning tree T of G, there is an edge not in T whose fundamental cycle with respect to T is a balanced 
separator. Goodrich |Goo95j observed that such an edge can be found by looking for an edge-separator in 
the cotree T* of T. 

We modify this approach slightly: let T* be a spanning tree of the planar dual G* of G such that T* has 
maximum degree 3. Let T be the cotree of T*, so T is a spanning tree of G. Root T* at an arbitrary vertex 
of degree one or two. Let t; be a leafmost vertex of T* such that the descendants of v comprise more than 
3/4 of the weight. Let e be the edge of T* connecting v to whichever child has greater descendant weight. 

Lemma 1. The fundamental cycle of e with respect to T is a balanced simple cycle separator. 

Proof. Simple algebra shows each of the two trees of T* \ {e} comprises between 1/4 and 3/4 of the total 
weight. One of these trees consists of the faces enclosed by G, and the other consists of the faces not enclosed 
by G, where G is the fundamental cycle of e with respect to T. □ 

Miller |Mil86j proved the following theorem (actually, a stronger version). 
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Theorem (Miller jMil86) ). For a planar triangulated biconnected graph G with weights such that the weight 
of each face, edge, and vertex comprises at most two-thirds of the total weight, there is a simple cycle C 
of length at most 2^j2\V{G)\ such that the total weight strictly enclosed is at most two-thirds and the total 
weight not enclosed is at most two-thirds. There is a linear-time algorithm to find such a cycle. 

In this paper, we do not use Miller's construction. We give another construction that, with the aid of some 
auxiliary data structures, can be carried out in sublinear time. For simplicity of presentation, we present a 
construction that achieves a balance of three-fourths instead of two-thirds, and refrain from optimizing the 
various constants that arise in our construction. 

2.5 Divisions 

Frederickson |Fre87j introduced the notion of r-divisions. Let G be an n-vertex planar embedded graph. 

Definition 3. A region R of G is an edge-induced subgraph of G. A boundary vertex of R is a vertex that 
is incident to an edge of R and incident to an edge of G not in R. 

Definition 4. An r-division of G is a collection of 0{n/r) regions such that each edge is in at least one 
region, and each region has at most r vertices and 0{y'r) boundary vertices. 

Frederickson's definition did not address the number of holes since it was not relevant in his algorithms 
(and in some subsequent algorithms building on his). The following definition follows the lines of that of 
Cabello |Cabl2j but our terminology is slightly different. 

Definition 5. A natural face of R is a face of R that is also a face of G. A hole of R is a face of R that is 
not natural. 

Definition 6. An i — division with few holes is an r-division in which every region has 0(1) holes. 

This differs from Cabello's definition in that his requires only that the average number of holes per region 
be 0(1). We use a stronger requirement because some algorithms depend on it. 
Our main theorem is as follows. 

Theorem 1. For a constant s, there is a linear-time algorithm that, for any biconnected triangulated planar 
embedded graph G and any r > s, outputs an r-division of G with few holes. 

2.6 Recursive Divisions and Decomposition Trees 

Some algorithms, e.g. the shortest-path algorithm of Henzinger et al. |HKRS97] . require that the graph be 
decomposed into regions which are in turn decomposed into regions, and so on. That algorithm requires 
roughly log* n levels of decomposition, so it would take more than linear time to find all the different divisions 
independently. We describe a simple decomposition of a planar graph that allows one to obtain such recursive 
divisions in linear time. 

Definition 7. A decomposition tree for G is a rooted tree in which each leaf is assigned a region of G such 
that each edge of G is represented in some region. For each node x of the decomposition tree, the region R^ 
corresponding to x is the subgraph of G that is the union of the regions assigned to descendants of v. 

Definition 8. A decomposition tree T admits an r-division with a few holes if there is a set S of nodes 
of T whose corresponding regions form an r-division of G with few holes. 

Theorem 2. For a constant s, there is a linear-time algorithm that, for any biconnected triangulated planar 
embedded graph G, outputs a binary decomposition tree T for G that admits an r-division of G with few 
holes for every r > s. 

Definition 9. For an exponentially increasing sequence r — (ri,r2, . . .) of numbers, a recursive r-division 
of G with few holes is a decomposition tree for G in which, for i — 1,2, . . ., the nodes at height i correspond 
to regions that form an ri-division of G with few holes. 

Theorem 3. There is a linear-time algorithm that, given a decomposition tree satisfying the condition of 
Theorem^ and given an increasing sequence r, returns a recursive r-division of G with few holes. 
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2.7 A Recurrence Relation 



The following lemma, whose proof is provided in the appendix, is used to analyze several quantities through- 
out the paper. 

Lemma 2. Let < P < 1, let c be a constant, and let 

T I \ ^ i + max{Q^} Tr{a.tn) if n > r 
^'■^"^ - \ ifn<r ' 

where the maximization is over {ai}f^i such that 
1. tti < 3/4 + c/ ^/n for i ~ 1, . . . , 8, and 

There exists a constant s such that for every r > s, Tr{n) is 0{n/r^-'^). 

2.8 Dynamic Trees and Euler-tour Trees 

The dynamic-tree data structure of Sleator and Tarjan jST83| represents rooted forests so as to support 
topological operations Link and Cut (removing and adding edges to the forests) and to support assignments 
of weights to nodes with operations that add a number to the weights of all ancestors of a given node, search 
among the ancestors for a node with minimum or maximum weight. Each operation can be performed in 
O(logn) amortized time. With care, given an initial tree, a dynamic-tree representation can be constructed 
so that the amortization is over 0(n/logn) operations. 

Many extensions and variants have been developed |ST83[ lEIT+92[ IFre97[ lABH+04[ I AHdLTOSi FTWOS) . 
For example, using, e.g., the self-adjusting top trees of Tarjan and Werneck |TW05j . one can search among 
descendants instead of ancestors, and one can represent embedded trees (adapting an idea of Eppstein et 
al. jEIT+92] ). „_ 

The Euler-tour-tree data structure of Henzinger and King [HK99| represents an Euler tour of a tree using 
a balanced binary search tree. It supports topological Link and Cut operations. In particular, for a vertex 
w of a tree, any consecutive (in the order used by the Euler tour) set of children of v can be cut in a single 
operation that requires a constant number of binary search tree operations, and hence takes O(logn) time. 
This is useful for representing embedded trees, by choosing the order in which nodes are visited by the Euler 
tour to be their cyclic order in the embedding. The Euler-tour tree can be decorated with additional labels 
to support ancestor and descendant searches in logarithmic time as well. 



3 Computing a Decomposition Tree 

We give a high-level description of the algorithm of Theorem [2] for computing a decomposition tree. The 
input graph G is assumed to be biconnected and triangulated. It follows that, for every region R of G, every 
natural face is a triangle. 

The algorithm is as follows. Given the input graph G, the algorithm performs some preprocessing 
necessary for the sublinear-time cycle-separator algorithm ( Algorithm [2|, and calls RecursiveDivide(G, 0). 
The procedure is given in Algorithm [l] 

This procedure, given a connected region R with more than s edges and given a recursion-depth param- 
eter i, first triangulates each hole of R by adding an artificial vertex and attaching it via artificial edges to 



each occurrence of a vertex on the boundary of the hole. Let R' be the resulting graph. See Figure 1(a) 1(b) 
for an illustration. The vertices and edges that are not artificial are natural. Triangulating in this way 
establishes biconnectivity of R' . 

Lemma 3. i?' is biconnected. 

Proof. Let u be a cut vertex of R. Let uv and uw be two consecutive (in the embedding) edges incident 
to u that belong to different biconnected components. Since vu, uw is the only v-to-w path in R, the edge 
vw is not in R. Hence uv and uvu belong to a face h of R that is not a triangle, hence h is a hole. After 
triangulation, vxh,XhW is another vertex-disjoint path from v to w, using the artificial node Xh- □ 
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Algorithm 1: RECURSiVEDiviDE(i?, i') 



1 let n = \V{R)\ 

2 if n < s then return a decomposition tree consisting of a leaf assigned R 

3 if £ mod 3 = then separator chosen below to balance number of vertices 

4 else if £ mod 3 = 1 then separator chosen to balance number of boundary vertices 

5 else if £ mod 3 = 2 then separator chosen to balance number of holes 

6 let i?' be the graph obtained from R as follows: triangulate each hole by placing an artificial vertex in 
the face and connecting it via artificial edges to all occurrences of vertices on the boundary of the hole 

r find a balanced simple-cycle separator C in R' with at most C\/n natural vertices 

8 let Fq, Fi be the sets of natural faces of R enclosed and not enclosed by C, respectively 

9 for i £ {0, 1} do 

10 let Ri be the region consisting of the edges of faces in Fi and edges of C that are in R 

11 Ti <-RECURSIVEDlVIDE(i?j, £ + 1) 

12 end 

13 return the decomposition tree T consisting of a root with left subtree To and right subtree 71 



Next, the procedure uses the SimpleCycleSeparator procedure (Section 4.4) to find a simple-cycle 
separator C consisting of at most c-^/n natural vertices, where c is a constant and n := is the number 

of vertices of R (which is equivalent to the number of natural vertices of R'). Depending on the current 
recursion depth £, the cycle separates i? in a balanced way with respect to either vertices, boundary vertices, 
or holes. 

Note that SimpleCycleSeparator is called on R' , which has more than n vertices since it also has 
some artificial vertices. However, we show in Lemma |6] that there are at most twelve artificial vertices, so 
even if we use a generic algorithm for finding a simple cycle separator, the bound of c^/n still holds for 
some choice of c (since n > s). In fact, our procedure SimpleCycleSeparator takes into account which 
vertices are artificial, and returns a separator consisting of at most 4-\/3n natural vertices, and possibly also 
some artificial vertices. The number of artificial vertices on the separator does not matter to the analysis of 
RecursiveDivide. 

The cycle C determines a bipartition of the faces of the triangulated graph i?', which in turn induces 
a bipartition {Fo,Fi) of the natural faces of R. For z = 0, 1, let Ri be the region consisting of the edges 
bounding the faces in Fi, together with the edges of C that are in R (i.e. omitting the artificial edges added 
to triangulate the artificial faces). See Figure l(c)[ 1(d) for an illustration. 



Lemma 4. // R is connected then Rq is connected and Ri is connected. 

The procedure calls itself recursively on Rq and Ri, obtaining decomposition trees To and Ti, respectively. 
The procedure creates a new decomposition tree T by creating a new root corresponding to the region R 
and assigning as its children the roots of To and Ti . 



3.1 Number of Holes 

The triangulation step (Line|6]) divides each hole h into a collection of triangle faces. We say a hole h is fully 
enclosed by C if all these triangle faces are enclosed by C in i?'. 

Lemma 5. Suppose that there are k holes that are fully enclosed by C. Then Rq has k + 1 holes. 

Proof. We give an algorithmic proof. See Figure [l] for an illustration. Initialize Rq to be the graph obtained 
from R' by deleting all edges not enclosed by C. Then C is the boundary of the infinite face of Rq . Consider 
in turn each hole h oi R such that a nonempty proper subset of /I's triangle faces are enclosed by C. For 
each such face h, C includes the artificial vertex Xh placed in h, along with two incident edges uxh and xuv 
where u and v are distinct vertices on the boundary of h. Deleting all the remaining artificial edges of h 
modifies the boundary of the infinite face by replacing uxh x^v with a subsequence of the edges forming 
the boundary of h. In particular, deleting these artificial edges does not create any new faces. 
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(a) A schematic diagram of a region R with four 
holes (white faces). 




(c) The region Rq consisting of the edges bounding 
the faces not enclosed by C together with the edges 
of C that belong to R. Equivalently, Rq is the 
subgraph of R' not strictly enclosed by C without 
any artificial edges and vertices. Rq has three holes. 



(b) The graph R' and a cycle separator C (solid 
red). Artificial triangulation edges are dashed (tri- 
angulation edges are not shown for the unbounded 
hole to avoid clutter). 




(d) The region Ri consisting of the edges bounding 
the faces enclosed by C together with the edges of C 
that belong to R. Equivalently, Ri is the subgraph 
of R' enclosed by C without any artificial edges and 
vertices. Ri has two holes. Note that a hole is not 
necessarily a simple face. 



Figure 1: Illustration of triangulating a hole and separating along a cycle. 



7 



Finally, for each hole h that is fully enclosed by C, delete the artificial edges of h, turning h into a face 
of R'q. The resulting graph is Rq, whose holes are: the holes of R that were fully enclosed by C, together 
with the infinite face of i?o- □ 

If the recursion depth mod 3 is 2, Line [t] of RecursiveDivide must select a simple cycle in R' that is 
balanced with respect to the number of holes. To achieve this, for each hole h of R, the algorithm assigns 
weight 1 to one of the triangles resulting from triangulating h in Line |6j then the algorithm finds a cycle C 
that is balanced with respect to these face- weights. 

Lemma 6. For any region created by RecursiveDivide, the number of holes is at most twelve. 
Proof. By induction on the recursion depth £, 

(A) if £ mod 3 = then R has at most ten holes; 

(B) if £ mod 3 = 1 then R has at most eleven holes; 

(C) if £ mod 3 = 2 then R has at most twelve holes. 

For £ — 0, there are no holes. Assume (C) holds for £. Since the cycle separator C encloses at most three- 
fourths of the weight, it fully encloses at most nine holes. By Lemma [sj Rq has at most ten holes. The 
symmetric argument applies to Ri. Thus (A) holds for £ + 1. 

Similarly, by Lemmajs] if (A) holds for £ then (B) holds for £+1, and if (B) holds for £ then (C) holds 
for ^ -1-1. □ 



3.2 Number of Vertices and Boundary Vertices 

If the recursion depth mod 3 is 0, Line [t] of RecursiveDivide selects a simple cycle in R' that is balanced 
with respect to the number of natural vertices. To achieve this, for each natural vertex v, the algorithm 
selects an adjacent face in R' , dedicated to carry w's weight. The weight of each face is defined to be the 
number of vertices for which that face was selected. Since each face in R' is a triangle, every weight is an 
integer between and 3. A cycle C is then chosen that is balanced with respect to these face-weights. 

If the recursion depth mod 3 is 1, the cycle must be balanced with respect to the number of boundary 
vertices. For each boundary vertex, the algorithm selects an incident face; the algorithm then proceeds as 
above. 

In either case, the total weight enclosed by the cycle C is an upper bound on the number of vertices 
(natural or boundary) strictly enclosed by C. Thus at most three-fourths of the vertices (natural or boundary) 
of R' are strictly enclosed by C in i?'. Similarly, at most three- fourths of the vertices are not enclosed by C 
in R'. 

The vertices of Rq are the natural vertices of R' enclosed by C (including the natural vertices on C, 
which number at most c^y\V{R)\), and the vertices of i?i are the natural vertices of i?' not strictly enclosed 
by C. Let n = \V(R)\ and, for i = 0, 1, let = \V{R,)\. We obtain 

no + ni < n + c\fn. (1) 

Moreover, if the recursion depth mod 3 is 0, then 

3 

max{noj"-i}^ -^'^^-c■^Tl. (2) 

Similarly, let 6 be the number of boundary vertices of i?, and, for z = 0, 1, let 5,; be the number of boundary 
vertices of Ri. We obtain 

5o + fei < 6 + c^/n. (3) 
Moreover, if the recursion depth mod 3 is 1, then 

3 

max{6o,^i} < -b^c\/n. (4) 
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3.3 Admitting an r— division 

Let N be the number of vertices in the original input graph G. Consider the decomposition tree T of G 
produced by RecursiveDivide. Each node x corresponds to a region R^. We define n{x) = \V{Rx)\- In 
this section, we show that, for any given r > s, T admits an i — division of G. We adapt the two-phase 
analysis of Frederickson |Fre87| . In the first phase (Lemmajsj), we identify a set of 0{N/r) regions for which 
the average number of boundary vertices is 0{y/r). However, some of the individual regions in this set might 
have too many boundary vertices (since the number of vertices and boundary vertices do not necessarily 
decrease at the same rate). We show that each such region can be replaced with smaller regions so that 
every region has 0{^) boundary vertices, and the total number of regions remains 0{N/r) (Lemma [9|. 

For a node x oi T and a set S of descendants of x such that no node in S is an ancestor of any other, 
define L{x, S) — —n(x) + X^j^es '^(j/)- Roughly speaking, L(x, S) counts the number of new boundary nodes 
with multiplicities when replacing x by all regions in S. 

Lemma 7. If y e Si then L{x,Si) + L{y,S2) = L{x,Si U S2 - {y})- 

Proof. Each node z G S2 contributes n{z) to both sides of the equation. Since y is in Si, it contributes n{y) 
to L{x,Si) (the first term on the left). It contributes —n{y) to L(y, S'2) (the second term on the left), and 
nothing to the righ-hand side. Every other node y' e 5*1 contributes n{y') to both sides. □ 

If the children of x are Xq and Xi, Equation ([T]) implies 

L{x, {xq, xi}) < c^/n{x). (5) 

Fix r and let Sr be the set of nodes y of T such that ?/'s region has no more than r vertices but the region 
of y's parent has more than r vertices. Note that no node in Sr is an ancestor of any other. Let x be the 
root of T- 

Lemma 8 (Total Number of Boundary Vertices). There are constants s and 7, depending on c, such that, 
for any r > s, L{x, Sr) 

Proof. For each node x of T, let £{x) denote the depth of x in T. For a node x of T, let Sr{x) denote the set 
of descendants y such that y's region has no more than r vertices but the region of y's parent has more than 
r vertices. For an integer n> r, define Br{n) to be mBxl L (x, Sr{x)) : x £ T,n{x) < n,i{x) = mod 3}. 
The lemma is a consequence of the following recurrence, which is justified below. 

Br{n) < 7 cy/n + maxy Br{ain), (6) 

{a.} Y 

where the max is over ai, . . . , ag between and f + ^ such that 



<l + 7c/^/i2 (7) 



and 

8 

^a,>l. (8) 

Let a; be a node such that r < \n{x)\ < n. Let j/i, . . . ,yfc be the rootmost descendants y of x such that 
!"-(?/) I < r OT £{y) mod 3 = 0. Note that k < 8. Repeated application of Lemma [T] and Equation ^ yield 

L{x,{yi,...,yk})<{k-l)cVn, (9) 

which implies that 

k 

^n{y,) <n+{k-l)cy/n. (10) 



9 



For 1 < i < fc, define ai so that aiU — n{yi). For k < i < 8, define ai — 0. Equation (10) then impHes 
Equation ([t]). Since each vertex of Rx occurs in at least one Ry., Equation Q holds. Let xo,xi be the 
children of x. For each 1 < i < k, Ry. is a subgraph of R^g or Rxi, so a^n < |n + Ci/n, which implies that 
ttj < I + c/y/n. 

By definition of Br{-), L{yi, Srijji)) < Br{ain), so we have 

8 8 

<^i?,(a,n). (11) 

2 = 1 1 = 1 

By Lemma [TJ 

k 

L{x,Sr{x)) < L{x,{yi, . . . ,yk}) + ^L(y„Sr{yi)), 

2=1 

so Equation ( 11 ) and Equation ^ together justify the recurrence relation (|6|. By Lemma[2] this recurrence 
relation yields that Br{n) is 0{n/y/r). □ 

Lemma [S] implies that X^^gs "-(2^) since the regions in Sr are disjoint except for boundary 

vertices, of which there are at most 0{N/y/r). For each parent y of a node x € Sr, the corresponding 
region Ry has more than r vertices, so the number of such parents is 0{N/r), so 15*^1 is 0{N/r). Let c' be 
a constant to be determined. For a node x, let S'r{x) denote the set of rootmost descendants y of x (where 
cc is a descendant of itself) such that Ry has at most c'^/r boundary vertices. Let S'^ = U^es i^'ri^)}- 

Lemma 9. The regions {Ry : y e 5*^} form an r-division with a constant number of holes per region. 

Proof. It follows from the definition of S^. that each region in this set has at most r vertices and at most 
c'^/r boundary vertices. It follows by induction from Lemma |4] that each of these regions is connected. It 
follows from Lemma|6]that each region has at most twelve holes. It remains to show that \Sl.\ is 0{N/r). 
For a node x of T, let b(x) denote the number of boundary vertices of R^. Lemma [s] implies that 

^ bix) < (12) 



We claim that, for every node x € Sr 



|5;(.t)| <max{l,^-12} (13) 



cwr 



Summing over x €z Sr and using Equation (12) then proves that \Sr\ is 0{N/r). 

We set c' — 40c. Proof of Equation (13) is by induction. If b{x) < c'^/r then 5*^(2;) = {x}, so the 
claim holds. Assume therefore that b{x) > d s/r.. Let yi,. . . ,yk be the rootmost descendants y of x such 
that b{y) < or t{y) - i{x) = 3, ordered such that |S';(yi)| > \S'riy2)\ > ■ ■ > |S'r(2/8)|- Let q be the 

cardinality of {i : \S'r{yi)\ > 1}. 

Case 0: q = 0. In this case, |S';(x)| < 8, and - 12 > 40 - 12 > 8. 

Case 1: q = 1. For some ancestor y of yi that is a descendant of x, the separator chosen for Ry is 
balanced in terms of boundary vertices. It follows that b{yi) < |6(x) + 3c^/r. By the inductive hypothesis, 

I^;(yi)|<^-12, so 

\S'r{x)\ = \S'riy,)\+k-l 

< ^„i2 + 7 

Cy'r 

U(x) + 3cJ^ 

< ^ ' ^ 12 + 7 

Kx) . „ ^ 

< ^--^+3-12 + 7 
c^r 4cy/r 

< ^-12. 

cJr 
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because b{x) > c'-^r — iOc^/r. 

Case 2: q = 2. In this case, b{yi) + 6(2/2) < &(a;) + Gcy^. Using the inductive hypothesis on yi and y2, 
we have 

= \SUyi)\ + \S:.iy2)\ + k~2 

< ^_i2+%l-i2-6 

b(x) + ^cJr „ „ „ 

< — 12-12-6 

< ^-12. 

cy'r 

Case 3: q > 2. This case is similar to Case 2. □ 



4 A Simple- Cycle Separator Algorithm 

In this section we present our cycle separator algorithm. As a one-shot algorithm, the input is a simple 
biconnected graph G with m edges and face weights, such that no face weighs more than 3/4 the total 
weight and no face consists of more than 3 edges. The algorithm outputs a simple cycle C in G, such that 
neither the total weight strictly enclosed by C nor the total weight not enclosed by C exceeds 3/4 of the total 
weight. The length of C is guaranteed to be at most 4^|£'(G)|. Since G is a simple graph, this implies a 
bound of on the length of the cycle. Better constants can be achieved, at the cost of complicating 

the algorithm and the analysis, which we avoid for the sake of simplicity and ease of presentation. 

In a similar manner, we aim for 3/4-balance to simplify the presentation. A balance of 2/3 can be achieved 
without significant complications. Note that handling face- weights can be used to handle vertex-weights and 
edge-weights as well; simply assign the weight of every vertex and every edge to an arbitrary incident face. 
This works as long as the resulting face weights satisfy the requirement that no single face weighs more than 
3/4 of the total weight. 

The cycle separator algorithm consists of a preprocessing step, which runs in linear time and computes 
certain auxiliary data structures used by the main procedure, SimpleCycleSeparator, which performs 
the computation. These data structures can be represented so that SimpleCycleSeparator takes sub- 
linear time. One auxiliary data structure is a tree /C, which is called the component tree. The tree JC 
captures the structural connectivity of dual BFS components of G. The dual BFS components satisfy a 



certain disjointness property (see Lemma 11). The other auxiliary data structures are a spanning tree T 



of G and its cotree T* . The spanning tree T satisfies a certain monotonicity property (see Lemma 12 1. 
The disjointness and monotonicity properties guarantee that the length of the cycle separator output by 
SimpleCycleSeparator is 4y^\E{G)\. 

Our algorithm for constructing the decomposition tree T invokes SimpleCycleSeparator multiple 
times, on various regions of G. It computes the auxiliary data structures /C, T, and T* once, for the input 
graph G, and efficiently updates their representation when separating a region into two regions. However, 
the disjointness and monotonicity properties mentioned above are slightly weaker (see Invariant [I] and In- 
variant [2]); They only apply to natural edges (disjointness) and to natural vertices (monotonicity). The 
implication is that the number of natural vertices on the cycle separator produced by SimpleCycleSep- 
arator is bounded in terms of the (squared root of the) number of natural vertices in the input graph. 
However, the cycle may also consist of some artificial vertices. As we argued in Section [3] artificial nodes on 
G do not affect the analysis. 

4.1 Levels and Level Components 

We define levels with respect to an arbitrarily chosen face /oo, which we designate as the infinite face. 

Definition 10. The level £^ (f) of a face f is the minimum number of edges on a /oo -to- f -path in G* . We 
use Lf to denote the faces having level i, and we use L^^ denote the set of faces f having level at least i. 
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Definition 11. For an integer i > 0, a connected component of the subgraph of G* induced by Ly^ is 
called a level-i component, or, for unspecified i, a level component. We use IC^i to denote the set of level-i 
components. A level-i component K is said to have level i, and we denote its level by t'^{K). A non-root 
level component is a level component whose level is not zero. The set of vertices of G* (faces ofG) belonging 
to K is denoted F{K). 

Note that we use K (not K* ) to denote a level component even though it is a connected component of a 
subgraph of the planar dual. K should be thought of as a set of faces. Thus we can refer to it as a subgraph 
of G* or of G. In the former case K is the subgraph of G* induced by the faces in the set K. In the latter 
case K is the subgraph of G induced by the edges that belong to faces in the set K. 

Lemma 10. For any non-root level component, the subgraph of G* consisting of faces not in F{K) is 
connected. 

Corollary 1. For any non-root level component K, the edges of 5g-'{F{K)) form a simple cycle in the 
primal G. 

In view of Corollary [ij for any non-root level component if, we use X{K) to denote the simple cycle 
in the primal G consisting of the edges of 5g'{F{K)). We refer to X{K) as the bounding cycle of K since, 
when viewed as a subgraph of G, K is exactly the subgraph enclosed by X{K). 

Lemma 11. Let K and K' be two distinct components. X{K) and X{K') are edge-disjoint. 

Proof. Let i be the level of K. The edges of X{K) are edges of 5g*{F{K)). Therefore, as edges of G*, they 
have one endpoint in K and one in a level-(j — 1) face. If an edge of Sg* {F{K')) has an endpoint in K then 
K 7^ K' implies the level of K' is at least i + 1, so it cannot be an edge of X{K). If, on the other hand, an 
edge of 5g* {F{K')) has an endpoint in level i — 1 then K ^ K' implies the other endpoint is not in K. □ 

The following definition is illustrated in Figm-e[2j 

Definition 12. The component tree K. is the rooted tree whose nodes are the level components and in which 
K is an ancestor of K' if the faces of K include the faces of K' . 

The root of the component tree is the unique level-0 component consisting of all of G* . 

Definition 13. An edge f f of G* has level i if f has level {i — 1) and f has level i. We write i^iff) for 
the level of ff . We use Lf to denote the set of edges of level i. 

Note that not every edge of G* is assigned a level. 

Definition 14. Let Lj denote the set of vertices of the primal graph G that are endpoints in the primal 
graph G of edges in Lf . 

Note that a vertex of G can be an endpoint of two edges at different levels i and j, so Lf and Lj are 
not necessarily disjoint. 

Definition 15. The level £^ (v) of a primal vertex v is defined to be min/ £^ (f) over all incident faces f . 
Note that LY is not the set of vertices with level i. 

4.2 The Primal Tree 

The algorithm maintains a primal spanning tree T. We start by describing the initial value of T and its 
properties. 

Lemma 12. There exists a spanning tree T such that, for any vertex u, £^ {parentrp[u)) < £^ (u). T can be 
computed in linear time. 
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(a) A triangulated graph (blue vertices and edges) along with (b) The corresponding com- 
a dual BFS tree (in red, faces at different BFS levels are indi- ponent tree K,. Each node 
cated by different shapes). The (primal) components K>k a-re of the component tree corre- 
indicated as shaded subgraphs. The deeper the level of a com- sponds to a connected com- 
ponent, the darker its shade is. ponent of G. 



Figure 2: Illustration of the component tree /C. 
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Proof. For a primal vertex u with £^ (u) — i, let f he a level-j face to which u is incident (ties are broken 
arbitrarily, but consistently). Let v,w be the other two vertices of that face /. The parent of u in T, 
denoted by parent7^(u), is the vertex in {v^w} with the smaller level (again, breaking ties arbitrarily and 
consistently). 

Since £^{u) = i, u is not incident to a level-(i — 1) face. Hence VjW must be adjacent to a level-(i — 1) 
face /'. Therefore, £^ (parent j.(u)) < i - 1 < (u) = i. 

To complete the definition of T, we choose an arbitrary vertex r incident to foo to be the root of T by 
assigning it to be the parent of the two remaining vertices at level i — 0. For convenience we set the level of 
the root vertex to —1. □ 

4.3 The Preprocessing Step 

To compute a simple cycle separator one has to first compute the component tree JC, the spanning tree T, 
and its cotree T* . This is done by the preprocessing step (Algorithm [2]) , which runs in linear time. 

Algorithm 2: Preprocessing(G) 

1 choose an arbitrary face as foo 

2 compute face, edge, and vertex levels £^ {■) , £^ {■) , £^ {■) , respectively 

3 compute the component tree K. 

4 initialize T to be the tree T as defined in Lemma 

5 initialize T* to be the cotree of T 

6 return {IC,T,T*) 
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The efficient implementation of RecursiveDivide, which computes the decomposition tree T in linear 
time, recursively separates regions of G. It maintains the component tree /C, spanning tree T, and cotree T* 
of the currently handled region R throughout the recursive calls. This is described in detail in Section [5] T 
is initialized to be the tree T of G. RecursiveDivide maintains the following invariants: 

Invariant 1. Let K and K' he two distinct components. X{K) and X{K') do not share natural edges. 

Invariant 2. f {v) < t {u) for any two natural vertices u and v of R such that v is an ancestor of u in T. 



Since G has only natural vertices and edges, Lemma 11 and Lemma 12 show that the invariants initially 
hold for R = G. 

4.4 Computing a Simple Cycle Separator 

We first provide an intuitive description of the separator algorithm SimpleCycleSeparator (Algorithm|3|. 
Let m denote the number of natural edges in G. 



1. The algorithm computes a balanced fundamental cycle G (Lines [3j|4|. If C consists of fewer than 
natural edges, then C is a short balanced simple-cycle separator. 

2. Otherwise, the fact that G is long implies that it intersects components at many (more than 2y^) 
consecutive levels of the component tree K.. (we say that C intersects a component K ii C has at 
least one edge in K — X{K) and at least one in G — K). The algorithm performs a binary search 
procedure on the range [I, h], where h are the minimum and maximum level £^ (v) of a vertex w e C, 
respectively. At each step of the binary search, we identify 

• the median level iq = [{I + h)/2\ , 

• the highest-level component K- intersected by C, whose level i- is smaller than zq and whose 
bounding cycle X{K-) has few (at most y/m) natural edges, and 

• the smallest-level component intersected by C, whose level is at least io and whose bounding 
cycle X{K+) has few natural edges. 
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The monotonicity of the primal tree T imphes that the number of natural edges of C between levels 
i_ and z+ (that is, the number of edges of C in — K+) is at most 2y/rn (see Lemma 15 1. 



3. If W{G — K-) > 3W/A we continue the binary search on the range [I, Similarly, if W{K-^-) > 3W/4 
we continue the binary search on the range [i+ , h] . 

4. Otherwise, the graph G' induced by the edges of X{K^),X{K+), and the edges of C in K- — K+ is 
a biconnected planar graph with at most natural edges, none of whose faces weighs more than 
3W/4. See Figure [sj The algorithm uses a simple greedy procedure, GreedyCycleSeparator, to 
output a balanced simple cycle separator in G' . 



Algorithm 3: SimpleCycleSeparator(C/), G is the tuple {G,IC,T,T*) 

1 let m be the number of natural edges in G 

2 let W := W{G) 

3 compute e* to be the 3/4-balanced edge separator of T*; let e be the primal of e* 

4 let C denote the fundamental cycle defined by e ^ T and T 

5 if C has at most 4-y/TO natural edges then return C /* C short enough */ 

6 let I, h be the minimum and maximum level £^ (v) of a vertex v £ C 

7 while I < h do /* binary search for iq in rainge [l,h) */ 

8 set all zo, i-, i+ to / + [(/i — /)/2j 

9 repeat /* sequential search for level i_ < iq with small boundary */ 

10 I z_ := i_ — 1; let be the unique component at level i- that C intersects 

11 until io < I or X{K^) has at most ^/m natural edges 

12 if W{G \K^) > 3W/A then h := i^; continue 

13 repeat /* sequential search for level > with small boundary */ 

14 I + 1; let be the unique component of level 1+ that C intersects 

15 until > h or X{K^) has at most ^Jm natural edge 

16 if W{K^) > 3W/4: then I i+; continue 

17 let G" be the graph induced on G by X{K^) U X{K+) U (c n [K^ - K+)] 

18 return GreedyCycleSeparator(G') 

19 end 



Note: To avoid clutter, the pseudocode of SimpleCycleSeparator does not handle the boundary case 



where the loop in line 10 terminates without finding a level with |X(X_)| < ^/m. In this case there is 

is 



12 



no need to shortcut G at a small level. Specifically, G is entirely enclosed by K^, the condition in line 
considered false, and X{K^) is considered to be an empty set of edges. A similar statement applies to K+. 



Lemma 13. The components and Kj^ defined in Lines 10 and 14 are well defined. 



Proof. Since I < i- < h, the cycle G must intersect some component at level z„. The monotonicity of the 
tree T implies that once a rootward path leaves a component at some level it never enters any component 
at that level again. The lemma follows since G is comprised of two rootward paths plus one edge. □ 

Lemma 14. For I < i < h, let Ki be the unique component at level i that is intersected by G. If the weight 
not enclosed by is greater than 3W/A then there exists a level I < io < i- such that the weight not 
enclosed by any with I < i < i^ is at most iW/A, and such that the weight enclosed by any Ki with 
io < i < h is at most iW/A. 

Proof. Let Ki be the unique level-^ component that encloses G. Since G is a balanced separator, the weight 
enclosed by Ki is at least W/A. Hence the weight not enclosed by Ki is at most iW/A. Let K be the 
component with maximum level that is intersected by G and whose enclosed weight is at least W/A. Let io 
be the level of K. By the above argument I < io, and since the weight not enclosed by is greater than 
3W/4:, io < i-. By choice of io, any Ki with io < i < h encloses at most W/A weight. □ 
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Figure 3: Illustration of the cycle separator algorithm: the fundamental cycle C is shown in red (solid and 
dashed). The level cycles X{K^) and X{K+) are indicated. The component if+ is shaded yellow. G" is 
the graph consisting of the edges of C that belong to — (solid red) as well as the edges of the cycles 
X{K^) and X{K+). 



A symmetric lemma applies to the case where the weight enclosed by is greater than 3VF/4. These 
lemmas show that the binary search procedure terminates. 

Lemma 15. Let G he a biconnected plane graph G with m natural edges and face weights, such that no 
face weighs more than 3/4 the total weight and such that no face consists of more than 3 edges. The 
procedure SimpleCycleSeparator finds a 3/ A-balanced simple cycle separator in G with at most 4-y/m 
natural vertices. 

Proof. Since G has face size at most 3, the fundamental cycle C is a 3/4-balanced simple cycle separator. 
If it consists of fewer than 'i^/rrl natural edges it is returned in Line [4j Otherwise, the lemma follows from 
the correctness of GreedyCycleSeparator (Lemma 16 1, provided that we show that G' is a biconnected 
subgraph of G with 0{y/rn) natural edges, none of whose faces weighs more than 3VF/4. We assume that both 
K- and K+ exist. The case where only one of them exists is similar. Consider the cycles X{K^),X{K-^-) 
and C. Since C intersects both K_ and K^, G' is biconnected. See Figure [3] for an illustration. 

To establish the bound on the number of edges in G', note that, by choice of X{K^) and X{K^), they 

consist of fewer than y/rn natural edges each. It remains to bound C n {K^ — K+) , which consists of two 

paths in T between X{Kj^) and X{K^). We claim that i+ — i_ < ^m. To see this, observe that by definition 
of i_ and 1+, for every level i_ < i < i+, the bounding cycle of the unique level-i component intersected by C 
consists of more than ^Jm natural edges. Since bounding cycles do not share natural edges (Invariant [l]) and 
since there are m natural edges in G, it must be that — i_ + 1 < ^/m. Combining this with Invariant [2J 



we conclude that 



Cr\{K_-K+] 



consists of at most 2 



natural vertices. 



As for the face weights, the weight of a face /' of G' is the total weight of the faces of G that are enclosed 
in G by the cycle /'. The faces of G' that correspond to the exterior of X{K_) and to the interior of X{Kj^) 
have weight at most 3VF/4 by the conditions in Lines 12 and 16 All of the other faces of G' are either 



enclosed (in G) by C or not enclosed by C. Therefore, the weight of each of these faces is at most 
since G is a balanced separator. □ 

Since any simple plane graph has at most three times as many edges as vertices, and since triangulating 
and making biconnected any plane graph G requires 0(|y(G)|) edges, we can restate Lemma 15 in a more 
general form that does not distinguish natural and artificial vertices: 
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Corollary 2. Let G be a simple biconnected plane graph G with face weights, such that no face weighs more 
than 3/4 the total weight, and such that no face consists of more than 3 edges. There exists a constant c 
such that the procedure SimpleCycleSeparator finds a 3 / A-balanced simple cycle separator in G whose 
length is at most cy''\V{G)\. 

We next discuss the procedure GreedyCycleSeparator that finds a simple cycle in a biconnected 
graph with no face weight exceeding iW/A. 

Lemma 16. Let G be a biconnected planar graph with m edges and face weights summing to W such that no 
face weighs more than iW/A. There exists an 0{m) algorithm that finds a balanced simple-cycle separator 
in G. 

Proof. The proof is constructive. Since G is biconnected, every face is a simple cycle. If W{f) > W/A for any 
face /, then / is a balanced simple cycle separator. Otherwise, use the following greedy algorithm. Initialize 
the set S with a single arbitrary face of G. Repeatedly add a face / to S, maintaining the invariant that the 
edges of Sc (S) (the boundary of S) form a simple cycle, stopping when W{S) first exceeds W/A. At this 
point, the edges of Sc-iS) form a simple-cycle separator. It remains to show that, in each iteration, there 
exists a face to add to S while maintaining the invariant, and to show how such a face can be quickly found. 

We claim that, if the boundary of 5 is a simple cycle, there exists a face f ^ S such that the boundary 
of S* U {/} is a simple cycle. We prove the claim using an inductive argument (illustrated in Figure [4]). Let 
e be an edge of S{S). Let /' be the face not in S to which e belongs. Let S" be 5 U /'. If S' has a simple 
boundary, we are done. Otherwise, S{S') consists of a set C of at least two simple cycles. Consider the graph 
G' obtained from G by deleting the non-boundary edges of S'. Each cycle in C consists of at least one edge 
of 6{S) and one edge of /', and bounds a connected subgraph of G' . Let H be one such connected subgraph. 
Since G is biconnected, any two vertices in H are connected by two vertex-disjoint paths in G. These paths 
can be transformed into vertex-disjoint paths in H by rerouting along the boundary of H . Such rerouting is 
possible since the boundary of is a simple cycle. Hence H is biconnected. Since H contains an edge of S, 
H Li S is biconnected as well. Now repeat the argument with G replaced hy S Li H. Since S Li H is strictly 
smaller than G (it does not include at least one edge — the one in /' n C" for some cycle G' E C that is not 
the one bounding H), the inductive argument proves the claim. 

We describe how, in each iteration, the greedy algorithm finds a face / such that the boundary of 5U {/} 
is a simple cycle. The algorithm maintains, for each face /, the number 7(/) of consecutive maximal subpaths 
of (5(5') on /. (A subpath may consist of a single vertex of S{S).) Note that S Li f has simple boundary iff 
7(/) = 1. The algorithm also maintains a list of the faces / not enclosed by S that have 7(/) = 1 and at 
least one edge shared with S. When a face / is added to the set S, the algorithm updates 7(-) as follows. 
Let S' denote S Li f. Since G is biconnected, / is a simple cycle. Let P be the subpath of / that consists 
of the boundary vertices and edges of 5" that do not belong to S. Let /i, /2, . . . be the faces other than / 
incident to edges or vertices of P in clockwise order along /. Note that the /^'s are not necessarily distinct 
(this is illustrated in Figure [5]). Let Pi be the maximal subpath (possibly a single vertex) of P along which 
fi is incident to P. For each subpath Pi, let and ef be the edges preceding and following Pi in S{S'), 
respectively. If both e~ and ef belong to decrease j{fi) by one. If neither e~ nor ef belongs to fi, 
increase j{fi) by one. Otherwise, 7(/i) remains unchanged. Since the work done when adding a face is 
proportional to the sum of degrees of the vertices of P, and since each edge or vertex belongs to P only once, 
when it first enters S, the total running time of this greedy algorithm is 0{m). □ 



4.5 Efficient Implementation 

It is fairly straightforward to implement SimpleCycleSeparator in linear time. In this subsection we 
describe the data structures that support implementing SimpleCycleSeparator in sublinear time. Our 
preprocessing procedure initializes the tree T to be the tree T (Lemma 12 1, the cotree T* to be the spanning 
tree of G* , rooted at foo whose edges are those not in T. T is represented using an Euler tour tree [ HK99) . 
T* is represented using a dynamic tree. It also computes the component tree /C and all level cycles in G. The 
component tree can be represented by a parent list. The level cycles are represented by splay trees |ST85| . 

For every edge e in G, let P'^{e) be the path in JC that consists of all components K such that e has one 
endpoint in X{K) and another endpoint strictly enclosed by X{K). We maintain these paths by storing the 
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Figure 4: Illustration of the proof of the auxiUary claim in Lemma |16| The boundary of the set S is shown 
in solid thick black. The face /' is shown in solid orange. H is a subgraph bounded by one of the simple 
cycles that comprise the boundary of 5' = S" U /'. A pair of vertices in H is indicated, along with two 
vertex-disjoint paths between them in the entire graph G. These paths can be transformed into disjoint 
paths in H by rerouting each of them along part of the simple cycle bounding H. The subgraph S U H is 
biconnected and strictly smaller than G. 




Figure 5: Illustration of the algorithm in the proof of Lemma [16] The boundary of the set S is shown in 
solid thick black. The path P is solid blue. The list of faces fi is /i = a, f2 = /3, fs — a, f4 ^ 7, — S. Note 
that the subpaths P3 and P4 both consist of the single vertex u. Before / is added to S, 7(a) — 1. Since e^, 
the edge of ^(5") preceding Pi, belongs to a, and e^, the edge of (5(5") following Pi, does not belong to a, 
Pi does not contribute anything to 7(0;). Since both ej and do not belong to a, P3 contributes +1 to 
7(a). Hence, after / is added to S, 7(a) = 2. 
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first and last components of P'^(e) in an array entry R(e). Note that R can be generated in linear time in 
the preprocessing step as part of the computation of the component tree. 

To implement the algorithm efficiently we need to support the following operations: 

• Line [l] - Find the number of natural edges in G. This can be done by keeping track of the number 
^artificial of artificial triangulation edges (Lemma |17| in Section [5] shows that this takes total 0(n) time 
over all calls to SimpleCycleSeparator, so we may assume m-artificiai is given). Let F be the number 
of faces in G. F can be found in 0(log |£^(G')|) amortized time by a single dynamic tree operation on 
T* . Since G is triangulated, the number of natural edges is given by 3F/2 — martificiai- 

• Line [2]- Find the total weight of faces in G. This is done using a single dynamic tree operation that 
returns the total weight of nodes (faces) in T* 

• Line|3]- Find a balanced edge separator in T* . This is done using a dynamic tree operation that returns 
a leafmost edge whose subtree has at least a certain total amount of weight {W/A in our case). This 
operation takes 0(log |i?(G')|) amortized time. 

• Line |4] - Find out if a path in T consists of at most a given number k of natural edges (in this case 
k — A,^frn). This can be done by traversing the path, and takes 0(A:log |_E(G)|) time (the number of 
artificial edges on a path is bounded by the number of holes of G, which, by Lemma [6j is constant). 

• Line [6] - Find least common ancestor in T. This can be done using the Euler tour representation in 
0(log |i?(G)|) amortized time. 

• Line[TOj [14]- Find leafmost ancestor in T with level at most i_. This can be done using the Euler tour 
representation in 0(log |i?(G)|) amortized time. 

• Find greatest-level component to which a particular edge belongs. This information is stored in the 
array R. 

• Line [12] [T6j [17] - Query the weight enclosed by a short simple cycle. This can be done by summing 
up the weight of subtrees as follows. For each cycle edge e, if e G T* and the rootward endpoint of 
e in G* is not enclosed by G, add the total weight of the subtree of e in T* to the sum. If e G T* 
and the rootward endpoint of e in G* is enclosed by G, subtract the total weight of the subtree of e 
in T* from the sum. Since querying the total weight of a subtree can be done using a dynamic tree 
operation in 0(log |i?(G)|) amortized time, and since the length of G is 0{^m), computing the total 
weight enclosed by G can be done in ©(^/rolog |i?(G)|) time. 

This shows that, assuming that the auxiliary data structures are given in the desired representation, 
SimpleCycleSeparator can be implemented in 0[^\E(G)\\o% \E{G)\) time. 



5 Maintaining the Representation of Regions Efficiently 

In this section we describe how to maintain the auxiliary data structures that represent the regions and 
are required by SimpleCycleSeparator throughout the recursive calls to RecursiveDivide. These 
include maintaining the embedding of regions, the primal spanning tree T and its monotonicity invariant 
(Invariant [2]) , the cotree T*, the component tree, and the level cycles. 

Consider an iteration of RecursiveDivide in which a region R is partitioned into two regions i?o and 
i?i along a cycle separator G. Note that, due to the recursive nature of the algorithm, we may assume 
that first the representations for i?o are obtained and the algorithm is invoked recursively on i?o- After the 
recursive call for i?o is completed, all changes are undone, and the representations for i?i are obtained. To 
achieve linear time for computing an r-division, it is required that the update be done in sublinear time 
in the size of i?. Essentially, the representation can be updated by performing only local operations on the 
vertices and edges of the separator G. This results in sublinear running time since G consists of 0(-y/|i?|) 
edges, and since each operation in the data structures we use takes 0(log|i?|) time. Lemma [2] with r — s 



4.5 



and 



shows that if the cycle can be computed in 0(|i?| ) time (for some /3 < 1), as established in Section 
if the representation can be updated in 0(|i?|^) time, then the total running time of the algorithm is linear. 
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5.1 Maintaining the Embedding 



The embedding is represented by using a splay tree |ST85| for each vertex u of i? to store the edges incident 
to V in cychc order. Conceptually, Rq is obtained from R by merging all the faces enclosed by C into a 
single face. Merging the faces enclosed by C corresponds to deleting all edges of R strictly enclosed by C. 
Similarly, Ri is obtained from R by deleting all edges not enclosed by R. We do not perform these deletions 
explicitly since that would not be efficient. 

We describe the procedure for Rq. The procedure for Ri is similar. We start by "cutting off" the edges 
that do not belong to Rq. To do this we reconstruct the cyclic order of edges around the vertices of the 
cycle C. Consider a vertex v of C. Let 7r„ be the cyclic permutation of edges around v in R. There are 2 edges 
of C that are incident to R. These cycle edges partition tTu into 2 intervals. The edges in one interval belong 
to Rq and those in the other interval do not. Using 0(1) cut operations, each taking 0(log|i?|) amortized 
time, we cut off the edges that do not belong to Rq. Since the length of C is 0{^/\R\) (it consists of 0{y/\R\) 
natural nodes, and at most 12 artificial nodes), the total time to update all vertices of C is 0(-\/|i?| log \R\)- 
Note that there is no need to explicitly delete the vertices enclosed by C from the representation because 
they are no longer pointed to from any vertex in Rq. 

Observe that the definition of Rq in Line 10 of RecursiveDivide and the definition of Rq in Line [6] 



of the recursive call of RecursiveDivide on Rq are equivalent to deleting all the artificial triangulation 
edges incident to C, and then triangulating the new artificial face of Rq with a new artificial node and 
new artificial edges. This can be done using a constant number of link and cut operations per node of C, 
each requiring 0(log|i?|) amortized time. The number of triangulation edges is bounded by the number of 
boundary vertices of Rq. Note, however, that even though our algorithm separates according to the number 
of boundary vertices every third recursive call, there is no guarantee on the ratio between the number of 
vertices in Rq and the number of boundary vertices of Rq. Instead, we bound the total time required for 
creating all of the triangulation edges throughout the entire execution of the recursive algorithm. Recall 
that for a node x of T, b{x) denotes the number of boundary nodes of R^. 

Lemma 17. '^^fz-j-b{x)log\Rx\ — 0{\Rx\), where x is the root ofT (i.e., R^ is the input graph G). 

Proof. Consider a node x of T. Let Cx denote the cycle used to separate R^. Observe that a vertex v on 
Cx appears as a boundary vertex in the regions of exactly two descendants of x in each level H > i{x) of T, 
unless V appears on another cycle separator Cy for some descendant y of x. Since the depth of the subtree 
of T rooted at x is bounded by ci log \Rx\ for some constant ci (the size of regions decreases by a constant 
factor every 3 levels) we have: 

^6(a;) < ^2ci|C,|log|i?,| (14) 

xer xer 

Note that if v does appear on Cy for some descendant y of x, its additional occurrences as a boundary vertex 
are accounted for by the contribution of Cy to the sum. 

Let T{n) be the maximum of X^ajeT^*"! I C*!; | log | | over all n-vertex triangulated biconnected planar 
input graphs. Let k, yi,. . . yk, ai, ... as be as in Lemma[8j Let Tx,yi,...,yk denote the subtree of T induced 
by nodes that are descendants of x but not descendants of any yi. Note that the number of internal nodes 
in Tx.yi....,yk is fc — 1, which is at most 7, and that for any internal node z of Tx,yx,...,yk, \Rz\ 5; \Rx\, and 
< \Cx\<c^\lQ. Hence, 

T{n) < 2ci • 7c^/n\og^ n + max V r(Q:,n) (15) 
Lemma [2] shows that T{n) is 0{n). □ 



5.2 Representation of the Spanning Tree T and its Cotree T* 

Let T and T* denote the primal and dual spanning trees of R, respectively, at the time the cycle separator C 
of R is found. In what follows we describe how to convert T into the spanning tree Tq of Rq. The procedure 
for converting T into the spanning tree Ti of i?i is similar. 

Let h be the new artificial face of Rq The tree T is cut at each vertex v oi h so that the only edges 
adjacent to in T are those that belong to Rq (i.e., not strictly enclosed by h). This breaks T into a forest F 
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in 0{yJ\R\ log \R\) time. Let Fq be the set of trees in F that consist of edges of i?o- The cotree T* is cut at 
the edges of h. This disconnects T* into a forest F* in 0(-\/|i?| log \R\) time. Let Fq* be the set of trees in 
F* that consist of duals of edges of i?o. 

The algorithm connects the forest Fq into a spanning tree Tq of i?o while maintaining the monotonicity 
invariant Invariant [2]) . 

Lemma 18. Let D be a cycle of natural vertices in R. Let v he a vertex of D such that the parent of v inT 
is enclosed by D, where enclosure is defined with respect to the root of T . There exists a vertex u £ D with 
£^{u)<£^{v). 

Proof. Since the root of T is not enclosed by D, there must be a vertex u G on the path in T from v to 
the root of T. By Invariant [2) (u) < £^{v). □ 

Recall that, after deleting all the triangulation edges incident to h, h is retriangulated triangulated with a 
new artificial vertex Xh and new artificial edges vx^ for each vertex of h. Let Vh be a vertex which minimizes 
£^ {v) among all vertices of h. The algorithm attaches Xh as a child of Vh in T. It then considers every 
vertex v of h. If v does not have a parent in _Fo (that is, the parent of w in T was enclosed by h in R, so 
after cutting T, w is a root of a tree in the forest Fq) then the algorithmsets Xh to be the parent of v in _Fo 
by adding the artificial triangulation edge vxh to Fq. This operation is performed by linking the tree of Fq 
whose root is v to x^. li v does have a parent in Fq we add vXh to Fq . At the end of the process, each edge 
of Rq is either in Fq or in Fq , and Fq contains no cycles. Hence Fq and Fq are spanning trees of Rq and its 
dual, respectively. 

Lemma 19. After the update, the tree Fq satisfies Invariant^ 

Proof. The only natural nodes whose parent is set by the update procedure are those on the new hole h of 



Rq that do not have a parent after F was cut at the vertices of h . By Lemma 18 for each such vertex v, 
£^ (vh) < £^{v). Since the algorithm sets Xh as the parent of v, and Vh as the parent of Xh, any natural 
vertex of i?o that is a proper ancestor of v is an ancestor of Vh- Invariant [2] is maintained since the ancestors 
of Vh are not changed by the update procedure. 

The number of operations required to obtain the updated representation of the tree and cotree for Rq is 



proportional to the number of nodes of h. By Lemma 17 the total update cost over all recursive calls to 



RecursiveDivide is 0(|G|). Obtaining the representation of the tree and cotree for Ri is similar. □ 
5.3 Representation of Component Boundaries 

The component tree K. is computed at the preprocessing step and represented by parent pointers. We do not 
update the component tree when RecursiveDivide partitions the graph. We only maintain and update 
the level-cycles (the boundaries of components) . Note that at any given recursive call /C correctly represents 
the enclosure relations between components in the current subgraph. However, /C also contains components 
that no longer belong to the current subgraph. This does not pose a problem since we only use /C for parent 
queries. 

For every edge e in G, let P'^(e) be the path in JC that consists of all components K such that e has one 
endpoint in X{K) and another endpoint strictly enclosed by X{K). The algorithm maintains these paths 
by storing the first and last components of P'^{e) in an array entry R(e). Note that R is generated in linear 
time at the preprocessing step and is be used when working with subgraphs at different recursive levels. 

We now describe how we update the representation of the level cycles. By construction, the cycle separator 



C intersects at most ^/m level-cycles (at most a single cycle at each level between i_ and ^4., see Lemma 13) 
These are the only level-cycles that need to be updated. We describe the procedure for LIq. The one for Ri 
is similar. As we mentioned above, these updates are first done for Rq. When the recursive call for Rq is 
completed, all changes are undone and we perform the updates for 

We first identify the components that are bisected by C. This is done by traversing the cycle C. For 
each edge e of G we mark each component K in P'^(e) whose level is at least i_ by the endpoint of e that 
belongs to X{K). This can be done by starting at the last component of P^{e) (Since G does not contain 
edges of any component with level greater than the level of the last component of P^{e) is at most 
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and following parent pointers in the component tree K. until either the other end of P'^ (e) is encountered or 
until a component of level i- is encountered. 

After traversing C, we have identified, for each marked component two vertices u and v on X{K). We 
cut the cycle X{K) at u and at v. Since cycle levels are represented using splay trees, this takes 0(log \R\) 
amortized time. This results in two paths, let Xq{K) be the one that consists of edges of i?o- The endpoints 
u and V of Xq (K) are incident to the new hole h oi Rq. We mark u and v to indicate that u can be connected 
to V using at two artificial triangulation edges (we do not explicitly add those artificial edges to form a 
cycle since artificial edges might be replaced with new ones at subsequent recursive calls). This is the new 
(implicit) representation of the level cycle. 

We note that level cycles are no longer disjoint. They may share artificial triangulation edges, but no 
natural edges. Thus, Invariant [T] is maintained. 



5.4 Representation of Vertex Weights 

Our algorithm separates regions according to three possible balance criteria: number of natural vertices, 
number of boundary vertices, and number of holes. In this subsection we discuss how these quantities are 
represented and maintained. 

As we mention in Section [3. 2 [ and in the beginning of Scction[4j vertex weights are handled by transferring 
the weight to some adjacent face. A cycle that is balanced with respect to the face weights is balanced with 
respect to the original vertex weights as well. 

We define three types of face weights: natural_weight, boundary _weight and hole_weight such that 
balancing according to each weight type is equivalent to balancing according to the corresponding quantity. 
Initially, all face weights are initialized to zero. For each vertex v of G, the algorithm associate with v an 
arbitrary adjacent face /, and increases natural_weight(/) by one. It maintains this vertex-face association 
in an array F so that, for every vertex, the associated face can be queried in constant time. Since the input 
graph G is triangulated, each face has weight at most 3, so no face has carries more than 3/4 of the total 
natural_weight (recall that the cycle separator is not invoked on graphs with fewer than s vertices). 

We now describe how to update the weights when separating a region R along a cycle separator C into 
subregions Rq and Ri. We describe the procedure for Rq. The one for Ri is similar. The algorithm checks, 
for each (natural) boundary vertex v of i?o if F[w] is a face of i?o- This can be done in 0(log |i?o|) time using 
the cotree T* for i?o. If ^[v] is not a face of i?o, the algorithm chooses an arbitrary face / of i?o that is 
adjacent to w, sets F[u] to /, and increases natural_weight(/) by one. This guarantees that each vertex of 
i?o assigns unit weight to one associated face of i?o- Note also that for each face / of i?o, all vertices assigning 
weight to / are vertices of i?o- Hence, balancing _Ro by natural_weight is equivalent to balancing by the 
number of natural vertices of Rq. The time to perform the update is 0(|C| log |i?|), which, by Lemma 17 is 
0(|G|) over the entire execution of the algorithm. 

To update boundary_weight, the algorithm sets boundary _weight(/) to one for every new triangu- 
lation face created when triangulating Rq. Similarly, to update hole_weight, the algorithm sets to one 
hole_weight(/) f or o ne arbitrary new triangulation face of the new hole formed when cutting R along C. 



Again, by Lemma 17 this takes total 0( G ) time. 
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A Proof of Lemma [2] 

Proof. We show by induction on n that, for n > r, Tr{n) < 5^n/r^~^ — 771^, where 7 > 1 and 6 > 1 
are constants to be determined. Note that S'fn/r^^^ — •fn^ is non-negative whenever 6n/r^^^ > , or 
equivalently, whenever n > ■ 

Suppose n > r, and assume the lemma holds for n' < n. Fix values of ai,...,ag. We show that 
pn^ + ^-Tr{ain) < Sjn/r^^'^ ~ jn'^ . Let I ^ {i : ain > gi/ll^ij-^ }, and let q = 

Case 0: q — Q. In this case, aiU < gi/i^i-fi) for all 1 < i < 8. Since > 1, it must be that r < n < 8r. 
Hence Tr{n) < pn^ , and 

S-fn/r^->^ - -fn'^ > S-/r^ - {S - S^r^ (16) 

Setting 7 > 8/9 we get 

Tr (n) < pn^ 

< Sprf^ 

< ir^ (17) 



Requiring (5 > 9 and combining (16 1 and (17) proves the inductive step in this case. 
Case 1: g > 1. For every i ^ X, 

r n 

CtiTL < — ^ < 



^1/(1-/3) - ^1/(1-/3) ■ 

Hence, for every i ^ I, ai < ^^^^.g) . It follows that 

We claim that S can be chosen sufficiently large so that J^iex'^i = 1 + e for some e > 0. To see this, 
assume without loss of generality that ai > a2 > • ■ ■ > ctg, and observe that 



^1/(1-/3) 



> «?-«i + l-7T7<^ (19) 



Let s be sufficiently large so that 3/A + c/y/n < 1. Since 1/8 < ai < 3/A + c/^/n, — ai is strictly positive. 
Define e' > to be ^ ai. Note that, by concavity of — a, e' is monotonically increasing with s. 



Substituting into ( 19 1 we get 

E"f^l + ^'-^ (20) 
Define e to be e' — gi/d-n-, , and set d sufficiently large so that e > 0. We have 

^af>l + e (21) 

For each i ^ I, Tr{ain) < 0. Hence 

pn^ + Tr{ain) — pn^ + (^5^ain/r^^^ — "fafn^^ 

i i^X 

< pn^ + 6jn/r^^^ '^^^ aj — jn^ 

< pn^ + (l + c/Vri)<57nA^"'^ -7n'^(l + e) 

= 5jn/r^^^ — 7ri^ + pn*^ + S"/Cy/n/r'^~^ — e^n^ 

= 6jn/r^''^ - 7n^ + + (57cni/2-'5/ri-'' £7) 

< 6jn/r^-^ -jn^ +n'^ {p + S-/c/r^^^ - ej) (22) 
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Let s be sufRciently large so that for any r > s, 6c/r^ ^ < e/2. Note that e does depend on s, but it since 
it is monotonically increasing with s, there exists such sufficiently large s. 



Then the right-hand side of ( 22 ) is at most 



S-fn/r'^''^ - 7^'^ +n'^ {p- 57/2) 



Setting 7 > p/e completes the proof. 



□ 
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